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GENERAL  OBJECTIVES  OF  TASK  R 


Many  of  the  wide  variety  of  problems  associated  with  the  use  of 
materials  at  high  temperatures  occur  in  connection  with  advanced  propulsion 
systems,. and  of  thes^,  some  of  the  most  critical  and  complex  are  encountered 
in  high  performance  rocket  motors.  Present  trends  in  rocket  design  and 
propellant  formulation  resulting  in  gas  flows  of  increased  temperature, 
pressure,,  and  corrosiveness  may  be  expected  to  aggravate  the  materials 
situation.  It  seems  self-evident  that  future,  long-range  solutions  to 
these  problems  must  relay  on  more  sophisticated  approaches  and  brokder 
knowledge  of  the  behavior  of  materials  in  high  temperature  gas  environments 
than  is  characteristic  of  the  usual  *'quick-fix" ,  or  expensive  cut-^and-try 
test  procedures.  It  also  seems  likely  that  improvement  in  structural 
materials  themselves  (i.e.,  higher  strength  at  high  temperature,  higher 
melting  point,  etc.)  has  reached  the  point  of  diminishing  returns,  so  that 
various  other  subterfuges  must  be  tried. 

Research  performed  under  Task  R  at  the  Applied  Physics  Laboratory 
or  its  subcontracting  agencies  is  intended  to  provide  some  of  the  fundamental 
knowledge  necessary  for  a  rational  understanding  of  the  behavior  of  materials 
at  high  temperatures.  This  is  a  very  broad  and  very  complex  field  involving 
many  different  scientific  disciplines.  While  no  attempt  is  made  to  rigidly 
limit  the  scope  of  Task  R,  the  general  emphasis  is  on  appropriate  research 
in  the  flow  and  physical  chemistry  of  high  temperature  gases  such  as  are 
characteristic  of  advanced  solid  propellant  rocket  motors,  and  the  phenomena 
basic  to  heat  transfer  and  cooling  techniques  in  such  environments,  rather 
than  in  the  properties  of  materials  themselves. 


A.  A.  Westenberg 
Program  Coordinator 


summary 


Thermal  Conductivity  of  Gases  . 

High  temperature  gas  thermal  conductivity 
measurements  by  means  of  a  line  source  in  a  laminar 
flow  are  being  carried  out.  Some  improvements  in 
the  apparatus  are  described.  The  theory  of  the 
thermal  conductivity  of  binary  polyatomic  mixtures 
is  outlined.  Experimental  results  on  N2>  C02>  and 
N2  -  €©2  mixtures  are  given  over  a  wide  temperature 
range. 

Rocket  Nozzle  Fluid  Dynamics  . 

This  project  aims  at  careful  experimental 
measurement  of  the  flow  properties  in  a  rocket  nozzle 
employing  a  typical  solid  propellant.  Some 
additional  evidence  of  flow  separation  near  the 
throat  is  given,  in  that  a  nozzle  with  larger  radius 
of  curvature  of  the  throat  section  gave  a  smaller 
effect  than  the  unmodified  throat.  Results  of  pitot 
tube  boundary  layer  measurements  are  presented. 
Refinements  in  the  infrared  spectroscopic  equipment 
for  determining  the  nozzle  gas  composition  have  been 
continued  and  some  spectra  are  shown. 

Rocket  Nozzle  Chemical  Kinetics  . 

The  purpose  of  this  project  is  to  apply 
the  results  of  Rask  R  and  other  research  to  an 
analytical  study  of  complex  chemical  kinetics  in  the 
flow  of  real  propellant  gases  through  nozzles.  The 
equations  and  numerical  program  for  machine, 
calculation  are  presented  for  the  particular  case 
being  studied.  Some  of  the  difficulties  being  en¬ 
countered  are  mentioned,  and  the  experience  gained 
so  far  in  surmounting  them  is  described. 


I.  THERMAL  CONDUCTIVITY  OF  GASES 
(N.  deHaas  and  A.  A.  Westenberg) 

Objective 

The  analysis  of  the  problem  of  convective  heat  transfer  both  in 
the  laminar  and  turbulent  flow  regimes  requires  that  the  properties  of  the 
gases  involved  be  quantitatively  understood.  The  coefficient  of  thermal 
conductivity  is  one  such  property.  Measurement  of  thermal  conductivity 
of  gases  has  not  been  a  forsaken  field  of  research^  and  several  common 
gases  have  been  investigated  up  to  temperatures  approaching  reaction  motor 
temperatures.  The  primary  shortcomings  of  past  and  present  research  in 
the  study  of  thermal  conductivity  of  gases  for  application  to  such  heat 
transfer  problems  are;  a)  Inability  to  apply  present  experimental  methods 
to  temperatures  above  about  1000° K  (most  measurements  are  made  at  substan¬ 
tially  lower  temperatures  than  this),  b)  the  complexity  and  difficulty  in 
using  the  usual  apparatus  for  this  measurement  ( thermal -conductivity  cells) 
and  3)  the  lack  of  measurements  of  species  (and  mixtures)  prevalent  in 
reaction  motors.  Other  shortcomings  could  be  listed. 

Recent  measurements  of  molecular  diffusion  coefficients  -  a 
closely  related  transport  property  -  of  common  flame  gases  at  temperatures 
in  excess  of  1000°K  applied  a  simple  concept  (Ref.  1)  that  could,  in 
principle,  be  used  for  measuring  other  gas  transport  properties  (conduc¬ 
tivity  and  viscosity).  In  this  method,  a  point  source  of  a  trace  gas  is 
located  in  the  center  of  a  uniform,  heated  laminar,  jet  of  a  second  gas. 

A  series  of  gas  samples  removed  immediately  downstream  of  the  source  and 
quantitatively  analyzed  is  used  to  determine  the  molecular  diffusion  coef¬ 
ficient.  The  technique  has  good  precision  (+1-2%)  and  is  quite  adaptable 


to  high  temperatures  and  various  gas  types. 
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Adaptation  of  this  method  to  the  measurement  of  thermal  conduc¬ 
tivity  of  gases  is  an  obvious  extension.  The  source  of  trace  gas  is 
replaced  with  a  line  source  of  heat  and  the  temperature  rise  downstream 
of  the  heat  source  is  measured  with  an  appropriate  technique.  The  method 
is  simple,  should  have  good  precision,  and  can  be  used  with  a  variety  of 
gases  at  high  ambient  temperature  (approaching  flame  temperatures).  Thus, 
the  objective  of  this  research  is  to  develop  the  line-source  technique  for 
the  measurement  of  gas  thermal  conductivity,  and  to  use  the  method  of 
obtaining  data  over  a  wide  range  of  temperature  on  pure  gases  and  mixtures 
of  propulsive  interest. 

Apparatus 

Previous  reports  (Refs.  2,3)  in  this  series  have  mentioned 
the  use  of  an  absolute  low-speed  anemometer  to  measure  the  velocity  of 
the  gas  stream  in  the  region  of  the  thermal  wake  measurements.  For  sim¬ 
plicity,  it  was  desirable  to  use  the  same  wire  as  the  A.  C.  heat  source 
required  in  the  velocity  measurement  (Ref.  4)  as  well  as  the  D. C.  line 
source  for  the  thermal  wake  width  measurement.  For  various  reasons,  an 
A.  C.  frequency  of  100  cps  was  the  most  convenient  and  satisfactory  for 
pulsing  the  wire  during  the  anemometer  measurement.  This  meant  that  the 
wire  had  to  be  less  than  1  rail  diameter  ^n  order  to  follow  the  100  cps 
signal  adequately.  A  common  source  (Pt  or  Pt-107o  Rh)  of  0.5  mil  diameter 
was  used  satisfactorily  for  a  number>.of.iruns  rin;'the.  temperature- range'5294®K 
to  1050°K.  There  was  a  difficulty  however.  Although  provision  was  made 
for  tightening  the  wire,  a  sufficient  amount  of  tension  could  not  be  placed 
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on  the  1.5  in.  long  source  to  prevent  it  from  oscillating  at  the  higher 
temperatures.  Attempts  to  place  sufficient  tension  on  the  source  usually 
broke  the  wire,  and  considerable  t£me  was  lost  in  replacement. 

Consequently,  it  was  decided  to  use  a  stronger  1-mil  wire  for 
the  line  source,  and  a  separate  0.5  mil  source  for  the  velocity  determi¬ 
nation,  since  the  latter  would  not  have  to  be  more  than  a  few  millimeters 
long  and  therefore  would  have  much  less  tendency ;;tOi  oscillate;  ‘.The, 
velocity  source  was  made  from  Wollaston  wire  mounted  on  a  set  of  wire 
prongs  which  may  be  located  at  various  points  along  a  diameter  across 
the  jet.  This  source  has  provided  good  temperature  signals  to  the  hot-wire 
pickup  at  jet  temperatures  as  high  as  1064“K  in  N2.  The  pickup  is  a  2mm 
dis'polvedrportion  of  0.0001  in.  Pt  Wollaston  wire  which  is  also  mounted 
on  a  set  of  wire  prongs.  These  prongs  are  attached  to  the  three-dimensional 
movement  on  which  the  thermocouple  is  mounted. 

One  of  the  following  sections  gives  some  results  on  the  measure¬ 
ment  of  the  thermal  conductivity  of  mixtures  of  N2  and  CO2.  The  mixtures 
were  prepared  by  passing  the  separate  gases  through  calibrated  critical 
orifices  and  then  into  a  common  chamber  before  entering  the  furnace.  The 
mixture  proportions  were  continuously  monitored  with  a  thermistor  thermal- 
conductivity  cell  which  was  calibrated  to  give  the  mole  fractions  of  the 
component  gases. 

The  temperature  of  the  jet  gases  is  measured  with  a  0.001  in.  Pt, 
Pt-107„  Rh  thermocouple.  With  regulated  gas  flow  through  the  furnace  and 
regulated  current  through  the  two  furnace  coils,  the  jet  temperature  will 
remain  constant  to  +1°K  for  a  period  of  several  hours.  This  degree  of 
stability  is  more  than  adequate  for  the  measurements  that  are  required. 
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Theory  of  the  Thermal  Conductivity  of  Binary  Polyatomic  Gas  Mixtures 

This  section  presents  a  summary  of  the  theory  of  the  thermal 
conductivity  of  a  mixture  of  two  polyatomic  gases  and  specifies  certain 
quantities  used  ini  the  next  section. 

The  thermal  conductivity  of  a  mixture  of  two  polyatomic  gases 
as  derived  by  Hirschfelder  (Ref.  5)  is  given  by 
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In  this  equation  X^  and  X^  are  the  thermal  conductivity  and  the  trans¬ 
lational  (i.e.,  with  internal  degrees  of  freedom  frozen)  thermal 
conductivity  of  the  mixture,  respectively,  X^  and  X?  are  similarly 
defined  for  the  pure  components  i  (i  =  1,2),  x^  are  the  mole  fractions, 
D,  „  is  the  binary  diffusion  coefficient,  and  D. .  are  the  self-diffusion 
coefficients  of  the  pure  components. 

The  translational  thermal  conductivity  of  a  mixture  has  been 
given  by  Muckenfuss  and  Curtiss  (Ref.  6;  see  also  Ref.  7)  and  is  as 
follows  for  two  components; 
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In  these  equations  are  the  molecular  weights,  T  is  the  temperature, 

•ff  ^ 

P  is  the  pressure,  and  4^2  2^2  usual  dimensionless  ratios  of 

collision  integrals  (Ref.  8)!;  for  the  particular  intermolecular  potential 

function  being  used.  The  Lennard-Jones  (6-12)  potential  was  assumed  here. 

In  the  calculation  of  X  for  mixtures  of  N.,  and  C0„,  X°  and  D.  . 

m  2  2'  i  ii 

were  computed  from  viscosity  data  on  the  pure  components  (Ref.  9)  using 
\°  -  (15/4)(R/M^)  (3) 

”11  -  <e4/5<>i)  1i 

* 

where  T)^  are  the  viscosities,  the  densities,  and  the  usual 
collision  integral  ratios  for  the  pure  components.  The  pure  component 
thermal  conductivities  were  the  experimental  values  obtained  in  this 
work,  while  the  binary  diffusion  coefficients  Dj^2  been  measured 
earlier  in  this  Laboratory  (Ref,  10).  The  ratios  of  the  collision 
integrals  were  obtained  for  the  Lennard-Jones  (6-12)  model  using 
(e/k)^2^=  (Ref..  10),  (€  /k)^p^  =  Z13°K,  and  (e  /k)^^  =  79.8°K  (Ref.  8). 

Results 

Measurements  of  the  thermal  conductivity  of  nitrogen  and  carbon 
dioxide  by  the  line  source  technique  in  the  temperature  range  from  300^ 
to  about  1000°K  are  given  in  Figs.  I-l  and  1-2.  Results  by  other  investi¬ 
gators  are  also  given.  It  will  be  seen  that  the  agreement  is  generally  good, 
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although  at  the  higher  temperatures  there  appears  to  be  an  appreciable 
discrepancy  for  CO2.  Whether  this  is  real  or  not  remains  to  be  seen  - 
the  data  are  being  checked. 

Fig.  1-3  shows  similar  data  for  a  50-50  mixture  of  N2  and  CO2. 
Mixtures  containing  25%  and  75%  N2  are  also  being  measured  and  the  data 
are  nearing  completion.  Fig.  1-4  gives  a  comparison  of  measured  and 
theoretical  values  for  the  binary  mixtures  over  the  range  of  composition 
at  two  different  temperatures.  The  theoretical  curves  were  obtained  as 
iri  the  previous  section,  while  the  experimental  values  were 
taken  from  the  vs.  T  curves  of  the  appropriate  mixtures.  It  is  ap¬ 
parent  that  theory  and  experiment  are  in  good  agreement  at  550“K  but  not 
at  room  temperature.  The  significance  of  this,  if  any^  awaits  further 
checking  and  study  before  any  interpretation  will  be  attempted. 
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II.  ROCKET  NOZZLE  FLUID  DYNAMICS 
(F.  K.  Hill  and  H.  J.  Unger) 


Oblectlve 

The  dynamics  of  high-speed  gas  flows  in  rocket  nozzles  has 

become  Increasingly  importar  :  to  the  efficient  utilization  of  high- 

energy  solid  propellants.  The  flow  characteristics  of  the  gas  mixture 

are  dependent  on  the  thermodynamic  properties  of  the  components  and 
certain  time-dependent  phenomena  such  as  chemical  reaction  rates, 

recombination  and  condensation.  Acquisition  of  the  experimental  data 

required  for  an  accurate  description  of  the  gas  flow  from  a  relatively 

simple  propellant  and  rocket  nozzle  combination  has  been  the  first 

objective. 

In  order  to  provide  realistic  and  representative  testing  con¬ 
ditions  for  the  experiment,  a  double-base  solid  propellant  (ARP)*  has 
been  chosen  which  provides  the  basic  constituents  comnon  to  most  high- 
impulse  propellants.  These  are  hydrogen,  nitrogen,  carbon  monoxide, 
carbon  dioxide  and  water  vapor.  Impurities  of  the  order  of  2  to  3%. are  • 

present  as  is  the  case  in  all  propellants  and  it  is  possible  to  add  in  ^ell- 
controlled  amounts  and  sizes  solid  particles  for  future  extensions  of  the 
experimental  studies.  To  begin  with,  it  is  believed  advisable  to  keep 
the  problems  as  uncomplicated  as  possible  while  still  retaining  the 
fundamental  aspects  of  the  phenomena  under  investigation.  Combustion 
pressure  and  temperature  of  the  grain  are  1100  psi,  nominal,  and  2500°K, 
respectively,  end  burning  is  employed  for  either  10-second  or  3.0-second 
operation.  These  conditions  are  sufficient  to  introduce  measurable 
effects  due  to  the  variable  gas  properties  as  the  gas  expands  through 
the  nozzle,  gas  non-equilibrium  conditions  in  the  expanded  flow 

*  The  propellant  is  manufactured  and  proof-tested  at  Allegany  Ballistics 
Laboratory. 
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and  heat  transfer  and  nozzle  divergence  angle  effects.  Associated 
phenomena,  such  as  erosion  and  deposition  of  solid  particles,  are  also 
present  for  observation  providing  additional  data. 

It  Is  anticipated  that  as  the  work  progresses  higher  specific 
Impulse  propellants  will  be  utilized  providing  temperature  In  the 
3000®K  to  4000®K  range,  and  that  investigations  will  be  undertaken  of 
the  effects  of  metal  additives  to  the  grain  and  some  properties  of 
materials  as  they  are  affected  by  the  rocket  gases. 

Flow  Separation  Near  the  Nozzle  Throat 

In  the  previous  progress  report  /(Ref.  1)  the  results  of  two 
tests  investigating  separation  in  the  nozzle  throat  region  were  described. 
The  third  test  in  this  series  was  subsequently  run  to  measure  the  pitot 
pressure  distribution  with  the  modified  nozzle  throat  geometry.  The 
pressures  obtained  indicated  that  a  shock  still  existed  near  the  axis  at 
this  station,  but  that  the  pressure  rise  in  the  shock- free  region  was 
smaller  than  occurred  with  the  unmodified  geometry.  This  substantiated 
the  previous  inference  that  the  increased  radius  of  the  throat  areas, 
providing  a  more  gradual  change  in';curvature  to  the  conical  surface, 
reduced  or  eliminated  the  separated  region.  The  static  wall  and  pitot 
pressures  were  constant  throughout  the  run  as  compared  to  decreasing 
values  in  unmodified  nozzle  tests.  It  was  concluded  that  the  shock  which 
remains  in  the  modified  nozzle  was  most  likely  due  to  a  transition  to 
turbulent  flow  which  was  further  substantiated  by  the  deposition  of 
solid  material  on  the  wall  surface  downstream  of  the  throat.  These 
measurements  thus  clearly  demonstrated  that  the  arbitrary  matching  of  a 
given  conical  expansion  region  to  a  throat  contour  generated  by  a 
circular  arc  may  lead  to  appreciable  deleterious  effects  in  the  usual 
nozzle  configuration.  For  optimum  performance  and  the  minimizing  of  heat 
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transfer,  this  area  should  be  contoured  so  as  to  provide  a  minimum 
discontinuity  in  the  second  derivative  of  the  profile. 

Auxiliary  Settling  Chamber 

An  auxiliary  settling  chamber  with  a  molybdenum  cup  for 
obstructing  particle  movement  was  designed  as  described  in  the  previous 
quarterly  report  (Ref.  1)  and  tested.  The  purpose  of  this  chamber  was 
to  permit  more  quiescent  gas  conditions  to  obtain  before  entering  the 
nozzle,  to  eliminate  the  radiant  light  from  the  propellant  flame  for 
optical  studies,  and  to  reduce  the  presence  of  solid  particles  which 
had  led  to  erosion  and  roughening  of  test  models  and  probes.  Several 
tests  were  run  with  the  above  configuration  and  it  was  found  that  the 
system  provided  cleaner  flow  as  evidenced  by  a  much  reduced  erosion  rate 
on  probes  and  models  and  that  there  was  a  relatively  small  supply 
pressure  drop.  Some  degradation  in  total  temperature  appeared,  however, 
amounting  to  about  SOO^K. 

Boundary  Layer  Measurements 

In  order  to  provide  the  fundamental  data  for  heat  transfer 
determination  to  the  surface  of  rocket  nozzles,  it  is  necessary  to  know 
the  detailed  boundary  layer  characteristics  under  realistic  flow  conditions. 
Experimental  measurements  of  the  boundary  layer  were  begun  in  the  rocket 
tunnel  with  the  probes  and  instrumentation  developed  earlier  in  the 
program.  The  measurements  were  made  at  a  station  where  the  Mach  number 
was  4.5  and  the  Reynolds  number,  based  on  boundary  layer  thickness,  was 
5.2  X  10^,  The  flow  in  the  boundary  layer  was  predominantly  turbulent 
and  of  interest  for  direct  comparison  with  turbulent  boundary  layer 
measurements  made  in  hypersonic  wind  tunnels.  • 

The  data  obtained  come  from  three  tests  in  which  eight  pitot 
probes  were  located  peripherally  around  the  section  at  varying  distances 
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from  the  wall.  For  the  probes  close  to  the  wall  the  tips  were  flattened 
by  swaging  and  filing  to  provide  an  overall  dimension  of  0.024  in.  with 
a  minimum  wall  thickness  of  0.003  in.  The  entire  probe  and  shank 
assemblies  were  fabricated  from  molybdenum  stock  to  eliminate  any 
junctions  or  seals  within  the  test  section  which  would  be  exposed  to 
the  rocket  gas.  The  results  of  the  measurements  are  shown  in  Figures  II-l 
and  II-2  which  give  typical  pitot  pressure  data  and  the  resulting  Mach 
number  distribution,  respectively.  These  profiles,  especially  that 
shown  in  Figure  II-2,  indicate  the  predominantly  turbulent  characteristic 
of  the  flow  from  about  0.05  in.  outward  to  the  free  stream  and  the  laminar 
sub-layer  at  the  wall.  Experiments  are  now  under  way  to  survey  the  total 
temperature  distribution  in  the  boundary  layer.  With  the  total  temperature 
data  and  the  Mach  number  distribution,  the  static  or  stream  temperature 
and  the  velocity  profiles  may  be  calculated.  This  will  then  provide  all 
the  needed  information  to  compute  both  the  heat  transfer  and  skin  friction 
coefficients  at  the  wall.  It  is  planned  to  make  surveys  at  several 
axial  locations  in  the  nozzle  and  as  far  upstream  as  the  environmental 
conditions  will  permit. 

Spectroscopic  Development 

The  development  of  an  infrared  spectrometer  for  the  purpose  of 
making  quantitative  determinations  of  the  composition  of  the  gas  in  the 
high-temperature  rocket  tunnel  during  operation  has  continued.  In  out 
previous  report  (Ref.  1)  evidence  was  given  of  mechanical  vibration 
troubles  in  the  optical  system  of  the  spectrometer.  The  most  likely 
offender  appeared  to  be  the  Nernst  glower  support  which  consisted  of 
sheet  micaK  The  positioning  of  the  image  of  the  source  on  the  entrance 
slit  of  the  spectrometer  was  very  critical  due  to  the  size  of  the  image 


PITOT  PRESSURE  (PSIA) 


DISTANCE  FROM  WALL  (lnch«>) 

Fig.  1 1-1  PITOT  PRESSURE  DISTRIBUTION  IN  THE  BOUNDARY  LAYER 

Measurements  were  made  on  the  Inner  nozzle  surface  at 
expansion  ratio  of  34,  free-stream  Hach  number  of  4.5. 
Data  shown  were  obtained  at  7  seconds  in  each  test. 


MACH  NUMBER 


0  0.1  0.2  0.3  0.4  0.5  0.6 

DISTANCE  FROM  WALL  (inches) 


Fig.  I 1-2  MACH  NUMBER  DISTRIBUTION  OBTAINED  FROM  DATA  GIVEN  IN  FIG.  1 

Nozzle  diameter  at  this  station  4  inches,  total  boundary  layer 
thickness  0.35  inches. 
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and  the  curvature  of  the  slit.  The  small  diameter  (1/16")  of  the 
available  Nernst  glowers  and  the  necessarily  complicated  collimating 
system  made  It  Impossible  to  avoid  some  amplitude  modulation  of  the 
light  source.  The  simplest  solution  to  this  difficulty  was  to  Increase 
the  width  of  the  source j  and  this  has  been  done  by  Incorporating  a  water- 
cooled  globar  which  is  x  5".  The  inconvenience  of  having  to  water 
cool  the  source  mountings  is  offset  by  the  mechanical  stability  of  the 
light  output. 

It  has  been  noticed  that  during  tests  with  the  rocket  tunnel, 
when  the  spectrometer  was  in  operation  but  not  in  the  test  cell,  the 
airborne  noise  produced  electrical  noise  in  the  recording  system.  By 
using  a  tape  recording  of  the  rocket  noise,  it  was  possible' to  subject 
the  electronic  system  to  a  sound  field  up  to  128  db  with  a  60  watt 
amplifier.  This  was  sufficient  to  simulate  the  noise  observed  during 
test  runs  and  the  offending  parts  were  found  to  be^  the  first  two  stages 
in  the  tuned  video  amplifier.  By  shock  mounting  the  tube  sockets, 
weighting  the  tube  shields  with  sheet  lead,  and  selecting  tubes,  it  was 
possible  to  eliminate  all  of  the  airborne  noise  effect  at  this  sound 
level  which  was  much  higher  than  that  encountered  during  test  runs. 

Another  improvement  that  has  been  incorporated  in  the  spectro¬ 
meter  is  the  replacement  of  the  Au-Ge  detector  with  a  liquid  nitrogen-  ■ 
cooled  PbSe  cell  built  to  our  specifications  by  Santa  Barbara  Research 
Center.  This  detector  is  about  10  times  more  sensitive  at  4.2  n  and 
5  times  at  6.0|a  than  the  Au-Ge  cell  previously  used.  However,  if  one 
wanted  to  extend  the  wavelength  range  beyond  6  n  ,  the  Au-Ge  cell  is 
superior  to  the  PbSe  cell  out  to  9  h. . 

In  order  to  make  absorption  measurements  at  different  wave¬ 
lengths  with  equal  accuracy,  it  is  necessary  to  have  the  product  of  the 
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transmission  factor  of  the  spectrometer,  the  spectral  sensitivity  of 
the  detector,  and  the  spectral  energy  distribution  times  the  effective 
slit  width  at  the  detector  constant.  This  condition  does  not  i%ive  any 
of  the  conventional  conditions  such  as  constant  energy  at  the  detector, 
but  does  allow  one  to  make  quantitative  comparisons  at  the  same  wave¬ 
length.  Our  intention  is  to  compare  the  absorption  maxlnnim  of  a  band 
observed  across  the  tunnel  with  that  of  a  collected  sample  of  the  same 
gases  at  the  same  density  at  room  temperature,  as  well  asv:tirith  a.-.synthetic 
mixture. 

An  attempt  was  made  to  program  the  slits  by  cutting  a  cam  that 
would  give  constant  amplitude  over  the  entire  wavelength  range.  Figures 
II-3  and  II-4  show  the  results  of  this  programming  and  consist  of  the 
absorption  spectrograms  of  the  flushed  and  unflushed  spectrometer  path, 
respectively.  The  cam  was  designed  experimentally  by  using  a  circular 
disc  graduated  in  10®  steps  with  the  radius  of  the  disc  chosen  to  give 
the  desired  slit  width  at  4  u.  The  motor  drive  of  the  cam  system  was 
disconnected  so  that  the  cam  could  be  set  by  hand  at  any  desired  wave¬ 
length.  The  slits  were  then  opened  by  slipping  a  feeler  gauge  between 
the  disc  and  the  push  rod  that  opens  the  slits  simultaneously.  By 
plotting  the  feeler  gauge  readings  against  the  angular  position,  a  graph 
was  obtained  of  the  desired  cam.  The  first  attempt  to  make  the  cam 
failed  because  10°  steps  were  too  large  and  0.001"  accuracy  in  setting 
the  milling  machine  table  was  too  rough.  In  the  second  attempt  to  cut 
the  cam,  the  feeler  gauge  readings  were  taken  every  5°  to  the  nearest 
0.001"  as  before,  but  the  graph  was  smoothed  to  the  nearest  0.0001". 

The  milling  machine  table  motion  was  checked  V7ith  a  dial  gauge  in¬ 
dependently  of  the  drive  screw.  The  resulting  rough  cam  was  then 
smoothed  on  a  filing  machine  and  finished  with  10|i  diamond  dust  on  a 
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Fig,  I 1-4  AMBIENT  ATMOSPHERIC  ABSORPTION 
PbSe  detector  with  programmed  slits. 
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brass  lapping  tool.  The  final  polishing  was  done  by  hand  with  the  cam 
in  operating  position.  The  high  spots  were  rubbed  down  with  a  fine 
oil  stone  until  Fig.  II-3  was  obtained  and,  although  there  remain  a  few 
undulations  of  the  order  of  0.0001"  on  the  cam,  it  was  felt  that  the 
curve  was  sufficiently  flat  for  the  present  requirements.  The  Venetian 
blind  effect  on  the  spectrograms  is  due  to  the  wavelength  blanking 
signal  applied  to  the  C.R.O.  beam  every  10®  of  rotation  of  the  wavelength 
cam. 

Figure  II -5  is  a  picture  of  the  rapid  scanning  infrared 
spectrometer  set  in  test  position  with  rocket  test  section,  "D",  partially 
assembled  to  show  details  of  the  absorption  cell.  "A"  and  "B"  are  the 
massive  optical  benches  which  support  the  infrared  source  with  its 
collimating  unit,  "C",  and  the  spectrometer  components,  respectively. 

"E"  is  the  combustion  chamber  with  its  loading  breech,  "F",  removed.  "G" 
is  the  housing  for  half  of  the  collimating  system  which  focuses  an  image 
of  the  light  source  on  the  entrance  slit  of  the  Perkin-Elmer  monochromator, 

"H".  "I"  is  the  housing  for  the  chopper  disc  which  cuts  the  light  beam 

» 

just  in  front  of  the  entrance  slit  at  1980  cps.  "J"  is  a  shield  over  the 

gear  and  cam  assembly  that  oscillates  the  Littrow  mirror  and  programs  the 

slits.  "K"  is  the  detector  assembly  which  houses  the  elliptical  mirror 

that  focuses  the  image  of  the  exit  slit  on  the  detector.  "L"  is  the 

sliding  brass  tube  for  focusing  the  light  on  the  detector  housed  within. 

On  the  top  of  this  tube  is  a  plastic  funnel  for  filling  the  Dewar  with 

liquid  nitrogen.  "M"  is  the  transistorized  preamplifier  for  the  video 

signal.  All  other  amplifiers  and  recording  equipment  are  located  in 

another  building.  The  "lectrodrier" ,  "N",  is  used  to  dry  the  air  that 

flushes  the  entire  system.  This  removes  the  atmospheric  CO  and  H  0  from 

^  2 

the  light  path  except  in  the  absorption  cell.  In  actual  operation  the 


Fig.  II-5  THE  RAPID-SCANNING  INFRARED  SPECTROMETER  IN  TEST  POSITION 

ON  THE  OPTICAL  DENCHES 
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the  CaF*^  crystal  window  housings,  "0",  are  connected  to  the  collimating 
boxes  with  flexible  rubber  tubes. 

After  all  these  improvements  Test  No.  125  was  run  with  the 
axis  of  the  absorption  cell  at  station  3.5"  from  the  nozzle  exit.  All 
previous  tests  had  been  made  at  station  5.0"  with  0.8"  circular  holes 
through  the  cylindrical  wall.  At  this  station,  fogging  of  the  windows 
reduced  the  light  transmission  as  much  as  50%  during  a  run  and  the  side 
tubes  were  always  coated  with  soot.  In  the  new  test  section  the  openings 
in  the  conical  wall  at  station  3.5"  were  made  5/16"  wide  and  1/2"  high. 

This  was  done  to  reduce  the  possibility  of  having  strong  shockwaves 
across  the  absorption  cell  and  to  allow  for  better  flt^shlng  of  the  CaF2 
windows . 

Figure  II-6  is  a  spectrogram  taken  during  this  test  and, 
although  the  noise  persists,  it  has  changed  in  character.  During  the 
run  there  was  no  window  fogging  and  inspection  after  the  run  showed  no 
soot  in  the  tubes. 

The  window  fogging  problem  was  solved  but  the  noise  problem 
lingers  on.  Before  considering  the  noise  problem,  let  us  look  at  the 
absorption  spectrum  and  compare  it  with  previous  tests. 

In  Figure  II-6  the  absorption  of  the  4.2  l-iC02  band  is  essentially 
the  same  as  Figure  III-IO  of  Ref.  1,  but  the  4.6l-iC0  band  is  extremely 
weak.  No  explanation  of  this  discrepancy  is  offered  at  this  time. 

Figure  II-7  shows  an  absorption  spectrogram  immediately  after 
burn-out  while  still  flushing  the  CaF^  windows  with  dry  nitrogen.  The 
various  oxides  of  nitrogen  are  observed  due  to  the  decomposition  of 
the  cellulose  acetate  inhibitor  casing  of  the  propellant.  It  was  noted 
after  the  run  that  a  part  of  this  case  had  melted  away  after  burn-out 
and  that  copious  fumes  of  NO2  (a  red-brown  gas)  were  issuing  from  the 


Fig.  I I -6  THE  ABSORPTION  SPECTRUM  OF  EXHAUST  GAS  IN  THE  WIND 

DURING  TEST  NO.  125 

Absorption  cell  at  station  3. 5"  on  the  25°  conical  nozzle. 
T  =  69C  °K;  p  =  4.6  x  10“^  Ibs/ft^; 

=  2.5  Ibs/in^;  M  =  4.42. 
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Fig.  II_7  ABSORPTION  SPECTRUM  IMMEDIATELY  AFTER  BURN-OUT 
AT  ATMOSPHERIC  PRESSURE  IN  THE  TUNNEL 

The  nitrogen  oxides  are  from  the  disintegration  of  the 
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combustion  chamber. 

1 

In  Figure  II-6  the  two  branches  of  the  CO2  band  at  4.2|i  r 
overlap  due  to  temperature  broadening  toward  the  long  wavelength  side 
sc  that  they  appear  as  one  band.  The  broadening  of  the  band  Is  obvious 
when  comparing  it  with  the  one  in  Figure  II-4  taken  at  ambient  temperature. 
In  an  article  by  C.  Tingwaldt  (Ref.  2)  the  absorption  of  this  band  was 
investigated  at  temperatures  300®  -  1000*K  and  the  broadening  of  the 
band  was  plotted.  Measurements  taken  from  Figures  II-4  and  II-6  show 
that  the  band  broadens  from  4.40 p  at  300®K  to  4.52 p  which  corresponds 
to  700®K  according  to  Tingwaldt.  The  theoretical  free-stream  temperature 
at  the  observation  station  is  690®K. 

The  noise  problem  has  been  pursued  and  has  been  further 
isolated  to  the  slit  mechanism.  Hammering  the  optical  bench  produces 
highly  damped  oscillations  characteristic  of  opening  the  slits.  Tapping 
the  base  of  the  spectrometer  casting  near  the  slit  assembly  produces  the 
same  effect,  so  it  appears  that  the  vibration  problem  is  one  or  two  steps 
nearer  to  solution. 

No  noise  has  been  attributed  to  the  detector  and  its  mirror 
system.  There  was  some  apprehension  as  to  how  stable  this  detector  would 
be  in  this  severe  environment  since  it  has  an  F/3  radiation  shield.  The 

r 

manufacturer  did  not  guarantee  the  stability  of  the  shield,  but  the  latest 
version  (#3)  has  the  inner  thimble  of  the  Dewar,  on  which  the  detector 
and  shield  are  mounted,  supported  by  a  plug  of  fibre  glass  between  it  and 
the  outer  glass  wall.  The  Dewar  is  held  in  a  fibre  bakelite  tube  by  two 
thin  rubber  washers.  This  makes  a  remarkably  stable  unit.  It  is  felt  that 
the  vibration  problem  will  soon  be  solved  so  that  quantitative  measurements 
may  be  made. 
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III.  ROCKET  NOZZLE  CHEMICAL  KINETICS 
(A.  A.  Westenberg  and  S.  Favin) 


Oblective 

Aerodynamics  and  chemical  kinetics  are  strongly  coupled  in  the 
flow  of  most  propellant  gases  through  an  exhaust  nozzle.  Understanding  of 
this  interaction  and  the  ability  to  predict  its  effects  must  be  founded 
solidly  on  the  fundamental  chemical  processes  Involved.  In  all  cases  of 
practical  importance,  these  basic  processes  are  exceedingly  complex.  As 
a  result,  all  of  the  theoretical  work  which  has  been  done  in  this  field 
has  been  confined  to  more-or-less  idealized,  simplified  chemical  systems  - 
generally  of  the  diatomic  molecule  dissociation-recombination  variety.  Even 
in  these  cases  the  analytical  and  computational  problems  are  formidable,  and 
various  approximations  have  been  proposed  for  dealing  with  the  situation. 

The  utility  of  these  approximations  in  the  more  complex  flows  typical  of 
real  propellants  remains  to  be  established.  The  aim  of  this  work  is  to 
apply  the  results  of  other  phases  of  the  Task  R  program,  as  well  as  work  done 
elsewhere,  and  examine  the  complex  chemical  kinetics  occurring  in  the  flow  of 
real  propellant  gases  through  nozzles  with  a  view  to  determining  if  and  when 
various  simplifying  procedures  are  valid. 

Basic  Equations  for  Numerical  Calculations 

In  a  previous  quarterly  report  (Ref.  1),  the  equations  governing 
the  flow  in  a  rocket  nozzle  with  chemical  reaction  were  presented  in  their 
most  fundamental  form.  The  discussion  was  appropriate  to  the  particular 
nozzle  and  propellant  being  used  in  the  experimental  phase  of  this  program 
(see  Section  II  of  the  present  report).  It  was  shown  that  under  the  conditions 
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assumed  for  the  nozzle  inlet  (T  -  2500‘’K,  P  -  1150  psia)  ,  the  only  chemical 
species  likely  to  be  important  were  CO,  C02»  H2»  H2O,  H,  OH,  and  N2,  and 
that  these  would  participate  in  the  reactions 


I.  CO  +  OH 


CO,  +  H 


II.  H  +  H2O 


H2  +  OH 


-2 


III.  H  +  H  +  M  ^ _  H  +  M 


IV.  H  +  OH  +  M  ^ _  H  O  +  M 


The  rate  constants  for  forward  and  reverse  reactions  are  related  to  the 
equilibrium  constant  by  (for  example) 


k 


1 


where  and  k  ^  are  in  units  of  cm^  mole  ^  sec.  ,  k^  and  k_^  in  cm  mole 
sec  and  the  equilibrium  constants  are  defined  in  terms  of  partial  pressures 
(atm.)  by 


K 


1 


^H  ^  ^CO  ^OH’ 


K 

3 


The  other  constants  are  similarly  defined. 


o^lo. 
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The  basic  equations  appropriate  to  the  problem  are  as  follows 
(the  composition  variable  is  mole  fraction  divided  by  mean  molecular 
weight  of  the  mixture  M) : 
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(10) 


in 

A. 


P  V 


A. 


(11) 


These  eleven  equations  are  to  be  solved  for  the  eleven  dependent  variables 

(seven  P,  T.p^and  v)it  being  presumed  that  the  area  ratio  A/A^.  (t  ■  throat) 

is  a  Icnown  function  of  the  distance  coordinate  z.  The  mass  flow  per  unit  throat 

area  m/A^.  is  in  the  nature  of  an  eigenvalue  of  the  problem  as  we  shall  see. 

Its  value  must  be  assumed  for  each  trial  solution.  The  as  well  as  T  and 

P  are  assumed  known  at  the  inlet  (subscript^) . 

For  actual  numerical  solution,  the  equations  may  be  simplified 

further.  The  pressure  P  may  be  eliminated  from  the  problem  first  by 

differentiating  Eq.  (10)  with  respect  to  z  and  combining  with  Eq.  (8).  This 

e 

gives 


dv  =  -  R 

dZ 


^  r-  dF, 
T  )  .  1 

^  dz 


dT 

dz 


I  ^  J- 


(12) 


p  dz  _ 

If  the  indicated  differentiation  in  Eq.  (9)  is  carried  out  making  use  of  the 
fact  that  for  ideal  gases  molar  enthalpies  are  functions  of  T  only  so  that 


dH. 


dz 


dH. 

1. 

dt 


dT 

dz 


=  C. 
1 


dt 

dz 


where  C.  is- molar  heat  capacity  at  constant  pressure,  we  get 


dv  =  -  J 

dz 


dz 


Differentiating  Eq.  (11)  gives 
1 


A/A.. 

c 


d(A/Aj.) 


dz 


P 


+  L. 

dz  V 


dz 


?  =  0 
dz 


(13) 


(14) 
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The  velocity  term  v^v/dz^roay  now  be  eliminated: between  Eqs.((13)  arid ((12), 
using  (14),  with  the  result  that 


Then,  eliminating  v^v/dz^  between  Eqs.  (14)  and  (12),  we  get 


(15) 


2 

dz 

T 

a7a^ 

d  (A/Aj.) 
dz 

.  p  d(A/Aj.) 

Wi 

r  i 

1 

-  Tp  ^ 

2 

A/Aj.  dz 

i‘) 

i"  1 

c 

(16) 


The  working  equations  have  thus  been  reduced  to  nine  in  number  -  (1)  through' 

(7),  (15)  and  (16)  and  the  variables  P  and  v  eliminated.  The  sequence  of 

events  in  the  numerical  solution  may  be  summarized  as  follows:  (i)  At 

any  point  z^  in  the  nozzle  where  all  conditions  are  known,  Eqs.  (1),  (2),  and 

(3)  are  used  to  compute  the  derivatives  of  F_.,  F„  ,  and  F„  The  quantities 

CO  H2  H2U 

appearing  on  the  right-hand  side  of  these  equations  are  all  known  functions 
of  T  (rate  constants  and  equilibrium  constants),  z  (the  area  function  A/A^.)  , 
or  assumed  (the  trial  value  of  A^/m).  (ii)  From  these  derivatives,  new 
•  values  of  F  ,  F  „  at  the  next  point  z  +  1  are  obtained  by  a  suitable 
integration  scheme,  and  then  from  Eqs.  (4)  -  (7)  the  new  values  of  the 
remaining  F^.  (iii)  The  temperature  derivative  at  is  obtained  from  Eq.  (15) 
and  thus  the  value  of  T  at  z^  (iv)  Similarly,  the  value  of  dp/dz  at  z 


is  obtained  from  Eq.  (16)  using  the  dT/dz  from  Eq.  (15),  and  then  p  at  z 

n+1 

All  the  dependent  variables  are  now  known  at  the  next  point  and  the 
procedure  begins  again. 
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Throat  Criteridn 

As  mentioned  above,  the  quantity  A^/m  is  an  eigenvalue  of  the 
nozzle  problem,  i.e.  there  is  only  one  value  of  A^/m  which  will  allow  a 
smooth  supersonic  solution  for  a  given  set  of  inlet  conditions  (F^,  T,  and  P) , 
kinetics  constants,  and  geometry.  The  problem  then  arises  as  to  how  to 
correct  A^/m.  In  the  case  of  "frozen"  flow  (isentropic) ,  the  criterion 
would  be  that  the  Mach  number  at  the  throat  be  unity,  i.e. 

Mach  No.  =  _ y _  ;  ^  ^i  ^i 

IFj  -  R 

Similarly,  in  the  "equilibrium"  flow  extreme  (also  isentropic)  a  Mach  number; 
can  also  be  defined  which  involves  an  effecMve  Y'.  The  latter  includes  the 
effect  of  the  shifting  composition  and  is  much  more  complex-  for  a  case  where 
the  mean  molecular  weight  is  not  constant  it  is  prohibitively  complex  for 
practical  use.  Therefore,  it  is  clear  that  the  concept  of  a  Mach  number 
defined  in  this  usual  way  is  not  a  natural  or  convenient  one.  Instead,  if 
one  returns  to  the  fundamental  equations  from  which  the  Mach  number  was  derived, 
the  desired  criterion  is  easily  found.  Thus,  for  any  flow  whether  isentropic 
or  not ,  Eqs.  (14)  and  (8)  may  be  combined  to  give 


d(A/Aj.) 

A/A^ 


dP  1 1  -  _dP 
2  \  dP 

p  V  \  • 


For  any  A^/m  above  the  critical  value,  dP  =  0  at  the  throat  where  d  (A/A^.)  =  0. 
But  for  the  critical  A^/m  where  the  flow  passes  from  subsonic  to  supersonic  and 
dP  0,  the  relation. 


must  hold.  Thus  the  criterion  for  a  correctly  chosen  A^/m  is  one  for  which 
Eq.  (17)  holds  at  the  throat  (to  the  desired  precision).  It  will  be  recognized. 


of  course,  that  in  either  the  frozen  or  equilibrium  isentropic  flow  extremes 
the  right  side  of  Eq.  (17)  would  be  written  P/S  p)^  which  is  also 


-  24  - 


I 


identified  with  the  sound  velocity  under  these  conditions.  Eq.  (17)  is 
more  generally  useful,  however. 

Numerical  Solutions 

The  actual  example  of  a  complex  chemically  reacting  nozzle  flow 
being  studied  is  taken  from  the  solid  propellant  work  being  carried  out 
experimentally  in  another  phase  of  the  Task  R  program  (see  Section  II  of 
this  report) .  The  inlet  conditions  calculated  for  chemical  equilibrium  of 
the  given  propellant  combustion  products  at  a  chamber  pressure  of  1150  psia 
are 


P  =  1150  psia 

o 


0.016712684 


T  =  2500‘'K 

o 

P  =  0.0092947139  g  cm' 

o 

M  =  24.373784 


0.0050861457 

0.0094785010 

0.0049351280 

0.000014076600 

0.000040940000 

0.0047602110 


The  chemical  kinetic  rate  constants  used  were  (Ref.  1) . 

3  3-1-1 

kj^  =  10  exp  (-10000/RT)  cm  mole  sec 


IS  3  - 1  ■ 

k^  =  10  exp  (-25000/RT)  cm  mole  sec 


21  2 

k^  =  5  X  10  /T 


6-2  -1 
cm  mole  sec 


,  ,«23;„2  6  ,  -2  -1 

=  10  /T  cm  mole  sec 

The  corresponding  equilibrium  constants  were  computed  at  100" 

intervals  from  the  basic  thermodynamic  data  (Ref.  2). 

The  nozzle  geometry  was  somewhat  simplified  mathematically  from 

that  reported  previously  (Ref.  1)  but  remained  essentially  the  same  for  all 

practical  purposes.  The  contraction  section  was  taken  to  have  a  radius 

governed  by  an  equation  of  parabolic  form 

r  =  0.2591090  (z  -  2.702431)^  +  0.8001031 


(0  ^  z^  2.702431  cm) 

For  the  expansion  section,  conical  geometry  was  assumed  of  either  25"  or 
12.5"  total  angle 

r  (25°)  =  0.2216900  (z  -  2.702431)  +  0. 8001031 ] 


I  (z >2. 702431  cm) 


r  (12.5")  =  0.10952  (z  -  2.702431)  +  0.8001031 


The  area  ratio  function  A/A^  was  then  computed  from  these  radius  functions 
(r^,  =  0.8001031  cm).  The  mathematical  discontinuity  at  the  throat  caused  by 
the  two  separate  radius  functions  led  to  no  difficulty. 

The  equations  outlined  in  the  preceding'  discussion  were  programmed 
for  solution  On  an  IBM-7090  computer.  The  simultaneous  first  order 
differential  equations  were  integrated  by  an  available  routine  employing 
a  fourth  order  Adams-Moulton  method  (variable  step-size)  with  Runge-Kutta 
starting  procedure.  Unfortunately,  the  complete  non-equilibrium  problem  has 
so  far  resisted  actual  machine  solution,  the  reason  being  not  yet  clear. 

The  difficulty  is  in  getting  the  solution  started  at  the  inlet,  where,  after 
one  or  two  stepfe,  the  integration  always  "blows  up".  A'-great  deal'offeffort 
has  been  spent  in  trying  various  cures  for  this  malady  without  success  so 
far.  It  seems  likely  that  the  trouble  lies  in  the  fact  that  one  is  looking 
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initially  for  very  small  changes  in  the  variables  -  all  of  which  are 
strongly  coupled.  It  was  thought  at  one  time  that,  because  the  mixture  was 
started  from  an  equilibrium  condition  where  all  the  bracketed  quantities  . 
in  Eqs.  (1),  (2),  and  (3)  are  zero,  the  difficulty  was  caused  by  extreme 
sensitivity  near  this  singularity.  However,  attempted  solutions  where  the 
initial  conditions  were  arbitrarily  selected  gas  mixtures  not  at  equilibrium 
did  not  behave  any  better.  Furthermore,  simplified  partial  solutions  such 
as  will  be  described  can  be  started  from  a  state  of  chemical  equilibrium 
with  no  problems  at  all,  so  that  this  does  not  appear  to  be  the  basic 
difficulty.  Work  on  this  is  continuing. 

The  iteration  procedure  for  finding  the  eigenvalue  k^/m  may  be 
illustrated  by  the  results  for  the  frozen  flow  case.  Iteration  is  started 
at  high  values  of  the  parameter  A^/m,  which  represent  mass  flov;  rates  below 
critical  and  thus  all  give  valid  solutions  (subsonic  throughout,  of  course) 
with  smooth  approach  to  the  throat.  A  series  of  such  gradually  decreasing 
values  of  A^/m  gives  results  for  frozen  flow  as  shown  in  the  following  table: 


A^/m  (cm  sec  g  ) 

2  ^cm) 

v^VdP^dP 

1.773100  x'lO"^ 

2.702431 

0.98853627 

1.773000 

2.702431 

0.99169568 

1.772900 

2.702431 

0  99759115 

1.772894 

2.702431 

0.99857523 

1.772892 

2.702431 

0.99911441 

Exact  eigenvalue 
»  lies  between 

1.772890 

2.7015790 

0.99998782 

these  two 

The  throat  lies  at  z=  2.702431.  As  A^/m  is  decreased,  the  value  of  the 
critical  quantity  vVdP/d-'p  approaches  closer  to  unity.  The  final  line  of 
the  table  shows  how  the  solution  starts  to  blow  up  when  A^/m  is  below  the 
eigenvalue,  since  it  is  not  possible  to  reach  the  throat  at  all  before  violent 
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oscillations  begin.  In  this  case  the  value  A^/m  =  1.772891  x  lO"^  would  be 
accepted  as  the  correct  eigenvalue  for  our  purposes. 

Once  the  correct  A^/m  is  determined,  there  remains  the  question  of 
getting  the  solution  through  the  throat  into  the  "supersonic"  regime.  The 
sensitivity  of  such  calculations  near  the  throat  is  well  known,  and  sometimes 
presents  a  problem  in  numerical  work  of  this  general  nature.  The  procedure 
we  found  most  successful  was  to  reduce  the  precision  index  of  the  integration 
routine  (a  parameter  which  represents  essentially  the  extent  to  which  extrapolated 
and  corrected  values  of  a  variable  are  expected  to  agree)  after  the  eigenvalue 
was  found,  until  the  solution  would  proceed  smoothly tthDou^h  .fche  ithcoat  .with 
steadily  decreasing  T  and  p  as  desired.  In  this  way  the  throat  region  was 
"skipped  overV,  so  to  speak.  If  desired,  a  more  precise  solution  in  the 
expanding  section  of  the  nozzle  could  then  be  obtained  by  starting  another 
calculation  at  a  point  slightly  downstream  of  the  throat  with  higher  precision 
index. 

In  addition  to-  the  complete  solution  to  this  nozzle  problem,  i.e. 
including  all  four  of  the  reactions  I  -  IV  as  non-equilibrium  processes,  it 
is  of  interest  to  examine  the  behavior  of  various  partial  solutions.  The 
main  effect  on  the  nozzle  flow  is  determined  by  the  degree  to  which  reactions 
I  and  II  take  place  durinf,  passage  through  the  nozzle.  The  net  effect  of 
these  two  may  be  written  as  the  water  gas  reaction. 

CO  +  H^O - 

which  is  exothermic  by  about  9.8  k  cal/mole  at  298‘’K.  Since  the  above  four 
species  (plus  the  inert  N^)  make  up  the  overwhelming  majority  of  the  total 
gas,  they  largely  determine  the  chemical  energy  effect  even  though  the 
recombination  reactions  III  and  IV  are  much  more  energetic  on  a  molar  basis. 

The  radicals  H  and  OH  are  present  in  very  small  amounts,  so  that  their  effects 
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are  only  indirectly  felt  via  their  influence  on  reactions  I  and  II. 

One  possible  simplifying  assumption  is  that  the  recombination 
reactions  III  and  IV  are  always  equilibrated.  Under  this  condition,  we  have 
the  relations 


K, 


"h^oVoh 


or  =  (F, 


H20^> 


(K3/P 


RTF„  ) 


(19) 


which  must  be  satisfied  throughout  the  expansion  process, 
of  III  and  IV  also  implies  equilibration  of  II,  since 


Equilibration 


Thus  the  only  rate  equation  remaining  is  that  for  reaction  I,  i.e.  Eq.  (1). 
Using  Eqs.  (18)  and  (19)  in  (1)  gives 


dF 

CO 

^  3/2  . 

P  •/ 

^CO  ^H20 

/  ■'3  1 

"CO2 

fSl 

%-j 

dz  ' 

*'/r  1 

(RT)^ 

L  *^4 

1  \l 

is  j 

This  equation  plus  the  element  conservation  equations  (4)  -  (7)  (F  and  F 

H  OH 

can  be  set  equal  to  zero  in  these  equations  for  all  practical  purposes)  and 
the  usual  Eqs.  (15)  and  (16)  are  sufficient  to  solve  the  problem.  This 
partial  solution  was  numerically  obtained  with  no  trouble  about  starting 
whatsoever. 

Another  solution  was  obtained  under  the  assumption  of  complete 
shifting  equilibrium  throughout  the  flow.  This  was  done  by  means  of  a 
standard  program  available  for  computing  equilibrium  composition  and  temperature 
after  isentropic  expansion  of  a  given  mixture  to  any  assigned  pressure  (Ref.  3). 
These  equilibrium  mixtures  were  then  related  to  area  ratios  and  nozzle  axial 
positions  z  by  appropriate  use  of  Eqs.  (9),  (10),  (11)  and  the  geometry 


(20) 


functions.  This  was  a  complete  equilibrium  solution  In  that  atoms  and 
radicals  (H  and  OH)  were  include^  unlike  the  earlier  equilibrium  solution 
which  considered  only  thu  water  gas  reaction  (Ref.  4). 
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